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SUMMARY 

Two sets of expressions for the short period perturbations due to 
lunar and solar gravitational forces in the orbital elements of artificial 
satellites are derived. The first set of expressions applies to the case 
of small eccentricity while the second set is applicable to orbits of 
arbitrary eccentricity. The results are tabulated in a form suitable for 
programming on high speed computers. 
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SHORT PERIOD LUNAR AND SOLAR PERTURBATIONS 
FOR ARTIFICIAL SATELLITES 


INTRODUCTION 

The influence of lunar and solar gravitational forces on the motion of a 
satellite were recognized shortly after the launchings of the first satellites 
(References 1 and 2). This led to analytic formulations for the perturbations 
in the elements. By expanding the disturbing function in powers of eccentricity 
and in terms of the mean anomalies of the satellite and perturbing body, accu- 
rate representations for the long period motion were obtained in Reference 1 . 

In addition, the results of Reference 2 are useful for satellites of small 
eccentricities . 

In order to accurately describe the short-period perturbations of satellites 
in orbits of moderate to high as well as low eccentricity, the results given in 
References 1 and 2 are both modified and extended in this report. 

In obtaining the solution for the case of low to moderate eccentricity, the 
expansion of the disturbing function in terms of the mean anomaly of the satellite 
as adapted in Reference 2 is carried out. The solution for satellites of arbitrary 
eccentricity is derived separately. Following the suggestion in Reference 2, the 
mean anomaly of the satellite is replaced by its eccentric anomaly. This permits 
replacing the infinite series in terms of the eccentricity of the satellite which 
may converge slowly, by polynomial in closed form. In addition, the convenient 
tabular form for representing the expansions is also adapted. 

However, rather than tabulating the disturbing function as done in References 
1 and 2, in this report the determining function (see Reference 3), and its differ- 
ential coefficients are given. In this way, the need for setting up the variation 
equations and performing the integrations by the reader is avoided. The perturba- 
tions in the elements are found by forming linear combinations of the tabulated 
coefficients. 


THE DISTURBING FUNCTION 

From Reference 4 the disturbing function, R , for the gravitational forces 
acting on a satellite due to the presence of a third body is given by 
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where 


k = Gaussian constant [k 2 = 6.670 x 10" S (cm 3 /kg. sec 2 )] 

m' = the mass of the disturbing body 

m' = 7.35 x 20 22 kg. 
m' = 1.99 x 10 30 kg. 

a, a' = the semi-major axis of the satellite and the disturbing body 
respectively. 

a' = 60.266011 earth radii for the moon 
a' = 23438.524 earth radii for the sun 

r,r' = absolute value of radius vector from the center of the earth 
to the satellite and to the disturbing body respectively 

C = cosine of the angle between r and r 

C is given by 


C = 


xx + yy + zz 
_ 


( 2 ) 


By substituting the expression for C given by Equation (2) into the disturb- 
ing function, R can be written in the form 


R = 




4 



.'2 


- 3 z 


'2 


.'2 


3 _ x 2 -y 2 
4 r 2 


.'2 


- y 


< 2 


.'2 



x y 3 yz 
.'2 ' „2 




1 0 ' 
r z x 


5x‘ 


.'2 


.'3 


r 2 y 


- 5y z 2 r ' 2 y *- 5y' z' 2 


.'3 


2 


1 3r 2 z -5z 3 3r <2 z' - 5z* 3 J> x 3 - 3xy 2 x' 3 - 3x'y' 2 
4 r 3 T 73 + 1 T 3 T 73 


5 3x 2 y - y 3 3x' 2 y' - y’ 3 15 x 2 z -y 2 z x' 2 z'-y' 2 z' 

8 r 3 r *3 4 r 3 r ' 3 

Next, the quantities enclosed in the brackets are expressed in terms of the 
orbital elements of the disturbed and disturbing bodies by the formulas 

— = c 2 cos(f + co + 9) + y 2 cos(f + co - 9) 


15 xyz x y z 


.'3 


( 3 ) 


— = c 2 sin(f + co + 8) - y 2 sin(f + co - 9) 


2 

— = 2yc s in(f + w) 
r 


( 4 ) 


x 


cos (f ' + co') 


=(l-2y' 2 )sin(f' + co') 
r' 


$ 

- 2 y' c' sin(f' + co') 


where 


c = cos 2 1/2 where I equals the angle of inclination of the orbit 
plane of the satellite to equatorial plane of the earth 


q = Q _ fi' where Q and fi' are the longitude of the ascending 
node of the satellite and of the disturbing body respectively 

y = sin 2 1/2 
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f = true anomaly of the satellite 

co = argument of perigee of the satellite. 

and where the primed elements refer to the corresponding elements of the dis- 
turbing body. 

The disturbing function may then be written in the form 

A' , A cos(qf + q'f' + qco + q' co' + v 6) 

A' vq , A vq cos (qf + q' f ' + qaj + q' co' + v 9) (5) 

In Equation (5) A vq is a function of c and y only while A' Vq , is a function of 
c' and y'. In the first summation of Equation (5) q and q' take on the five sets 
of values (0,0), (2,-2), (2,0), (0,2), and (2,2). In the second summation q and q' 
take on the eight values (1,1), (1,-1), (1,3), (1,-3), (3,1), (3,-1), (3,3) and (3,-3). 

Next the development proceeds in terms of the angular variables that are 
dependent upon the mean anomaly of the satellite or of the disturbing body. In 
Equation (5) first hold q, q' , and v constant and let 

a - qf + qco + q ' co' + v9 (6) 

to obtain a term of the form 






cos(a + q' f ') 


( 7 ) 


The part enclosed by the brackets may be expanded through the use of 
Cayley’s Tables (Reference 5) or by formulas from elliptic motion in Reference 
4 into a series of the form 




B. , q , cos (a + i'V) 


(8) 
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In (8) i ' takes on both positive and negative integers and zero values. Also 
B/ q * is a power series in the eccentricity e' of the disturbing body, which 
has a value of approximately .055 for the moon and .017 for the sun. 

The procedure outlined above may be continued for particular values of i'. 
Each term of the series, (8), may then be written in the form 


k 2 


m 


i'p 


a p A' , A B. , , 

pq Vq i'q' 



cos (q f + t) 


where t - i' V + qa> + q'w + v 6 


(9) 


For satellites in orbits of small eccentricity an expansion of the quantity in 
brackets in (9) converges rapidly. However, when the eccentricity of the satel- 
lite's orbit is large it may be more convenient to follow the suggestion of Ref- 
erence 2 and expand this quantity in terms of the eccentric anomaly of the satel- 
lite. Therefore, two distinct approaches to the solution will be pursued. The first 
case will be that of small eccentricity with the quantity in brackets in (9) expressed 
in terms of a truncated infinite series in eccentricity and mean anomaly, while 
in the second case that of arbitrary eccentricity this quantity is expressed in 
terms of a finite series of eccentric anomaly with exact coefficients involving 
the eccentricity. 

While it is true that the solution for the case of arbitrary eccentricity is valid 
for small eccentricity also, the solution for the small eccentricity case is much 
simpler. Thus, if the satellite for which the short period lunar and solar per- 
turbations are to be computed is one of small eccentricity considerable computing 
time and labor may be avoided by using the formulas for that case rather than those 
for the case or arbitrary eccentricity. 


THE DISTURBING FUNCTION FOR LOW ECCENTRICITY 

By means of formulas in Reference 5, (9) may be expressed in the following form 


^ aPA vq' A v q B iV T C iq cos(U + i'4‘ + q 

«p i 1 


a) + q'oj' + vff) 


( 10 ) 
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where C is a function of the eccentricity of the satellite. The coefficients A vq , 

A' », B i » q >, and C. are given in Tables 1 through 13 for low eccentricity. Thus, 
for this case the mean anomalies of both the satellite and the disturbing body are 
used in the expansion of the disturbing function and later the determining function. 
However, if i = i' = 0, then that particular term is not a short period term and 
should be omitted. Since it would be inconvenient to list the terms for which i = 

0 and i' =/ 0 and i / 0, i' = 0 separately from the tables, an asterisk (*) is placed 
next to the entry i = 0 or i* = 0 to remind the reader to omit the term for which 
both i and i' are zero. 


THE DETERMINING FUNCTION AND PERTURBATIONS 
FOR LOW ECCENTRICITY 

Let 0 denote the first order part of the determining function (Reference 3), 
then the perturbations in the Delaunay Elements L, G, H, ■i, g, hare given by 


8Lri2. 


Sg = -M 
3G 


Sh = - 


where 


L = ix ' a 


1/2 „ 1/2 


% - mean anomaly 


G = L(1 - e-*) 


2 \ 1/2 


g = argument of perigee 


H - G cos I 


h = longitude of ascend- 
ing node 


In this report g is replaced by co , while h and fi are equivalent. Also 


I = 1 - 2T 2 


From Reference 3, we find that 0 must satisfy the partial differential equation 


dQ 3Q 

n irt + n W R p 
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where 

n = mean motion of the mean anomaly of the satellite 
n' = mean motion of the mean anomaly of the disturbing body 
n' = 13.2/day for the moon (approximately) 
n' = .986/day for the sun (approximately) 

R p = Periodic part of the disturbing function. 

After integration, Equation (12) consists of two groups which are of the form, 


0 


k 2 


lei: 


A' , A B. , , C. 

Vq Vq l'q' lq 


sin (i't + i'V + vd + qco + q'o/) 


in + i'n' 


(13) 


In group 1 (p, p' ) = (2,3) and in group 2, (p, p' ) = (3,4). 

A typical term of R p is of the form 

A cos (i'i + i'-C' + 0) 

and the perturbations then resulting from such a term are obtained from 


8L = 


BQ 



SH = 


BQ 

B0 


g g = - 

Sh = - 


BQ _ 2L BQ BQ 

BL /x Ba 6 l 3 Be 

BQ _ G BQ cos I BQ _ G ^ -2'y 2 ) 

BG ^2 Be + G B(cos I) e j ^2 Be 4Gy By 

BQ _ J_ BQ _ J_ BQ 

BH G B(cos I) 4Gy By 


(14) 


(15) 
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In order to obtain (BQ/Ba) one must make use of 




dn 3 n 

da 2 a 

Therefore, if the typical term in R p is A cos (i^ + i.'V+ <f>) where 


( 16 ) 


k 2 r 


a p A' , A B , , C. 

vq' vq l q iq 


then the perturbation may be obtained by applying the following formulas 


BB 


A’ , 

A B 

, , c 


a?' 1 

Bl" 

a'P' 

vq' 

Vq 

i q 

" ( 

in + i'n' ) 2 

BB 

k 2 m 

A' 

B , , 

3C. 

iq 


1 

Be 

a'p' 

vq' 
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Be 

in 

+ i'n' 

BB _ 

k 2 m 

A' , 

BA 

Vq 

B , , 

c 

1 

By 

a'P' 

rt Vq' 

By 

i'q' 

»q 

in + i'n' 


[(2p + 3) in + 2pi'n'] 


(17) 


where 


B = 


A 

in + i'n' 


ILLUSTRATIVE EXAMPLE OF THE USE OF THE TABLES 
FOR LOW ECCENTRICITY 

In order to illustrate the use of the tables, the contribution to the perturbation 
in -i resulting from the term defined by the indices q = 3, q* = 3, v - 2, i =5, 
and i = 2 is obtained. From the first of Equations (15) 


s-e = - — ^ ° 2 

M 3a e L 3 ’ 

and from the first two of Equations (17), where A represents a single coefficient 
in the series for Q 
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From Table 12 


_ k2m A' A B C 

3 a — A 23 A 23 5 3 S 3 
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Therefore, due to this single term 


u= 

V M Ba e L 3 ° e 


] cos (2t + 5-i* + 3 oj + 3cJ + 


2d). 
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LOW ECCENTRICITY 


TABLE 1 
q = 0, q' = 0 
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LOW ECCENTRICITY 
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LOW ECCENTRICITY 
TABLE 3 
q = 2, q' = 0 
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LOW ECCENTRICITY 


TABLE 4 
q = 0, q' = 2 
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LOW ECCENTRICITY 
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THE DISTURBING FUNCTION FOR ARBITRARY ECCENTRICITY 


By means of formulas in Reference 4, (9) may be expressed in the following 
form 


, iilil' a p A„ A' , B., . y 1 C. cos (i E + i ' V + qw + q' a? + v 6) (18) 

,p ' ’■'q *V / , iq v n ' 

i 

where C. is a function of ft and where 

L. , iq .. -- ■ .. . 

/3 = - (1 - = e - (19) 

e i + vr^ 

The quantity ft is thus of the same order as e , the eccentricity of the satellite's 
orbit. Further, in C. q the expansions in ft are in terms of closed polynomials 
rather than infinite series as would be the case if the expansion had been in 
terms of the mean anomaly of the satellite instead of the eccentric anomaly , E . 

Thus the expansion of R, or even (r / a)R which as shown below is useful in 
deriving the determining function, can be expressed for given values of q and q' 
in the form 

cos (iE+ i'V + v& +qaj+q'w') (20) 


The coefficients A' v , A Vq B., , and C iq are given in Tables 1 through 13 for 
arbitrary eccentricity. q q q 

The expansion of (r/a)R is used in finding the determining function which is 
now discussed. 

THE DETERMINING FUNCTION AND PERTURBATIONS 
FOR ARBITRARY ECCENTRICITY 



The determining function for the arbitrary eccentricity case is obtained in a 
slightly different fashion from that of low eccentricity. The development for this 
case is the same through Equation (12) however. The difference in the form of 
the determining function is due to the need of using eccentric anomaly instead 
of mean anomaly in order to obtain closed expressions involving eccentricity. 


Since in this section the eccentric anomaly of the satellite rather than the 
mean anomaly is used. Equation (12) takes the form 


|2. + m i^Q = j_ R 

BE a 2V na p 


(21) 


23 



where 


m 



m* +M' 



3/2 


and where m ffl and M' are the mass of the earth and disturbing body respectively. 

But it has been shown that a typical term in the expansion of (r/a)R p is of the 
form 


A cos (iE + i ' V + 4>) ■ 

Hence we may write Equation (15) in the form 


+ m = — co s(iE + i' £'+<£) + m , cos E (22) 

oE cK/ n o'i 

where nij = me . Equation (22) is a type considered by Hansen, Reference 6, 
whose solution may be represented in the form 

Q = Q 0 +m 1 Q 1 +m2Q 2 + ... (23) 

Hansen's method for finding the expansion represented by Equation (23) has been 
described in Reference 2, which we repeat for the convenience of the reader. 

We first neglect the terms of order and higher in Equations (22) and (23) 
to obtain the equation 


+ m _S2. cos (i E + i' + 4>) (24) 

3E ~dV n 


The solution of Equation (24) is given by 


Q = Asin(iE+i , ^ < + <£) 
0 n(i + i' m) 


( 25 ) 


which may be verified by noting that Equation (25) satisfied Equation (24) . 


24 


The next step is to neglect terms of order m \ and write 


Q* *.>"(!»* i'-t' **> tm Q 

n(i + i m) 1 1 

which is substituted into Equation (22) to obtain 


( 26 ) 


BQ X 3Q X 

BE +m W 


BQ 

u 


2-cos E 


and in general we may write 


BQj , BQ. BQj.i 


BE ' m 'd-f 34* 


cos E ’ 


(27) 


(28) 


Then, from Equations (23), (25), and (27) the solution for Q may be found to order 


m 4 


Q = . (Q^ + meQ; + m 2 e 2 Q') 

n(i + l m) u 


(29) 


where Q‘ 0 ,Q' and Q' 2 are trigonometric functions given in Table 14. 

Substituting the value of Q from Equation (29) into Equation (13) we find that 




SG = 2 + me 

Ba> 

?Qo 

w 


SH = 


+ m e 


bq; 

+ m2e2 ^ 

a A 

BE 

BE/ 

r n i + n' i ' 

bq; 

2 2 BQ 2^ 
■f m * a * - i 

1 A 

Bo; 

i m tr J 

Bw / 

n i + n' i ' 

b q ; 

bq; 

, 2 2 2 

) A 

B<9 

m e dd 

/ $ • i 

/ ni + n l 


( 30 ) 
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The differential coefficients with respect to E , u , and 9 are listed in Table 14. 
Similarly 


U = 


BQ 

BL 


2L 30 | G 2 BQ 
M Ba e L 3 Be 


BQ G BQ + cos I BQ G BQ _ (l-2y ) BQ (31) 

BG ~ ^2 Be G B(cos I) eL 2 Be 4Gy By 


sn = - 22 = _ I bq _ l bo 

BH G B(cosl) 4 Gy By 
We also need the formulas 


BE 

3e 


a . _ 

-- — s i n E 


d/3 _ (i + ay 

de 2(l-/3 2 ) 

dn _ _ 3 n 
da 2 a 


(32) 


For group 1 of Equation (5), the quantity, A, appearing in Equation (29) is given 
in Tables 1-13 for Arbitrary Eccentricity. It is of the form 


k 2 m' 


,'3 


3 2 ^vq' 


A vq B i'q' 


C. 

iq 


so that BQ/Ba is given by 


26 


[Qq +meQ; + m 2 e 2 Q'] 


BQ 

Ba 


k 2 m < a(7ni + 4n' i') 
2a' 3 (ni + n'i') 2 


A' 


Vq 


A B. , , 

Vq i q' 


c. 

iq 


Bme Q! Bm 2 e 2 Q' 


ni + n' i' L Ba 


Ba 


(33) 


where BmeQj/Ba and Bm 2 e 2 Q' /Ba are given in Table 14. Similarly, for group 2 
of Equation (5) we find 


BQ _3k 2 m'a 2 (3ni+2n'i*) 

Ba 2 * 4 / . i • i \ o 

a (n i + n i )•* 


A' 


vq' 


Vq 


B i'q' 


C.q [Q' + meQj m 2 e 2 Q' ] 


ni + n' i ' 


Bme Q' 
_ Ba 


Bm 2 e 2 Q' 
Ba 


(34) 


For group 1, BQ/Be is given by the formula 


BQ k 2 m' a 2 f(l + /3 2 ) 2 3C ; 

= A.. ' A B ; , — — — 


2(1 ~/3 2 ) 1 q B/3 


[Qq + m eQj + m 2 e 2 Q'] 


n i + n' i ' >- 


^ a 
BE r 


bq; „ 

— sinE+mQ! + me — — — sinE+2m 2 eQ' 

- BE r 


+ m 2 e 2 — sin E 
BE r 


(35) 


For group 2, the factor k 2 m'a 2 /a' 3 is replaced by k 2 m'a 3 /a' 4 . Finally for 
group 1, BQ/3 y is given by 


BQ_ = k 2 m' a 2 

By" a , 3 *V 


BA 


Vq 


By 


B. . < C : . 


iq 


[Q; + me Qj + m 2 e 2 Q' ] 


(36) 
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The corresponding derivative for <3Q/ ty for group 2 is again obtained by replacing 
the factor k 2 m'a 2 /a' 3 by k 2 m'a 3 /a' 4 in Equation (36). 


SPECIAL CASES 

Several interesting special cases arise which simplify the calculation of 
certain terms in the determining function, they are, (i.) i i 0, i' = 0, (ii.) i = 0 
i' ^ 0 and (iii.) i = i* = 0. 

(i.) The Case i / 0, i' = 0 

The differential equation for Q, Equation (21), then becomes 


_ JL £_ R 
3E n a p 

Hence Q is of the form 


(37) 


Q = ^aPA' v A vq B 0q , C iq sin(iE + ^) 
na' p 

So that Q = Q 0 with i' =0 and thus agrees with the result of the general case, 
(ii.) The Case i = 0, iV 0 

The differential equation for Q, Equation (21), for this case becomes 

n' — = R (38) 

W p 


Now, from Equation (9) we find Q is of the form 


Q = 


k 2 m' 


V A 'vq* 


A „ B. . . C ft s in(i* V + 4>) 

Vq i q Oq v 


F or this case we find that 

Q = Q 0 

with the subscript i set equal to zero. 


( 39 ) 
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(iii.) The Case i = i' = 0 

Let R s be the terms of R that are independent of both i and i ' obtained by 
averaging R over the mean anomalies of the satellite and the disturbing body. 
Then Equation (21) may be written as 


§ = - < R - R ,> 

3E na s 


(40) 


Integrating we find 


R 

Q=_iL (E-'t), 
n 


(41) 


since 


di - — dE and 
a 



— RdE = R . 
a s 


Thus, (r/a)R = R s for those terms in (r/a)R for which i = i' =0. Thus 


R 


.'3 




l/q 


B 


0q' ^ Oq 


cos (qoj + q'oj' + vQ) 


+ 


k 2 m' a 3 \ 


iV 


Vq B 0q' ^0q 


cos (qaj + q'oj' + v9') 


(42) 


where the summations are over q , q' and v . We have 
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R 


a 



8G =. 



( 43 ) 


SH ---- R 


n 


s 


where R s is the expression obtained when the cosine terms of R are changed to 
sine. The perturbations in -t, w, and fi are given by Equations (21). The expres- 
sions for 3Q/B a, BQ/9 e , and >Q/ By are 


±Q 

"d a 


z.g.- k ^ (E - 1 ) Y 

2 n a ' 3 ^ — 1 


A ; q < A v q v c o q cos (q 


+ q' u' + v 0) 


9 a 2 k 2 m' 
+ 2n a' 4 


(E 


-<>L 


A' # A v „ B , C 


v q 


v q 0 q ' 0 q 


cos (q co + q' cJ + v 9) 


(44) 


3Q _ k 2 m' a 2 
' r)o na' 3 


a r 
— sin E 

r 


Y, A<l V Avt i Bo< »' C ° q 


cos (qco + q' &>' + v 0) 


k 2 m' a 3 


a _ 
— sin E 
r 


A' Vq , A^ B 0q , C flq 


cos (q co + q' co' + v 0) 


| k 2 m'a 2 (1 + /i 2 ) 2 
2 n a' 3 1 -/-l 2 




E 


BC 


A' 




A B , 

Vq Oq 


Oq 


'd/3 


cos (qo; + q'a/ + v 6*) 


30 


+ 


cos (q CO + q'oj' + v 6) 


k 2 m' a 3 (1 + / 3 2 ) 2 
2na’ 4 i-P 2 




Oq' 



( 45 ) 


BQ k 2 m' a 2 
na' 3 



BA 


y q 


By 


B 0q' C 0q 


cos (q^ + q'aj* + to 8) 


k 2 m' a 3 


n a 


M 


(E 


v— i BA 

-oZ/v-sv 


vq 


B 


Oq ' ''Oq 


C. cos (qoj + q'o/ + v 9) 


(46) 


ILLUSTRATIVE EXAMPLE OF THE USE OF THE TABLES 
FOR ARBITRARY ECCENTRICITY 

In order to illustrate the use of the tables, the contribution to the perturbation 
in co resulting from the term defined by the indices q = 2, q' = -2, v - - 1 , i'=-2, 
and i = 3 is obtained. 

It is assumed that the "mean" elements a, e , i , co and fi are known. Then, 
for a given value of E may be determined from Kepler's equation by using the 
mean eccentricity. Then r /a is known also as are the elements of the disturbing 
body. 

1. From Table 14 we find 


Q' = sin (3 E - 2 't 1 + 4>) where cp = 2 co - 2 co' -6 



sin (4E - 2l‘ + <t>) 1 

2 - m 2 


sin (2E - 2 l* + 0) 
1 - m 


O' - s * n (E - 2 V + 4 s ) sin (3E - 2 V + <j>) 
" 2 (1 -m) (1 - 2m) + 2(l-m)(2-m) 


sin (5E -2V +4>) 
2 (2 -m) (5 -2m) 



= 3 cos (3E - 2 V + 4>) 
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c)Q ^ 1 

— *-=-=-£. cos (4E _ 2 V + 0) _-±— cos (2 E - 2 +0) 

3 E 2 - m 1 - m 


B Q 2 cos (E-2 V +<t>) _ cos (3E-21* + </>) 5 cos (5E - 2 V + <fi) 

BE 2(1- 2m) (1 -m) (1 - 2m) (2 - m) 2 (2-m)(5-2m) 


2. From Table 2 we find 


Then 


Also, 


A:j_ 2 - 3c' y' 


B,-, 


A-i2 = cy ; 


C 32 = 


(1+/3 2 ) 3 


. _ k 2 m' a 2 . . 

71 ^~ i -2 ^-12 ®3-2 ^3 


dC 32 - 1 + 5/3 2 


d/3 (1+/3 2 ) 4 


We then find dQ/de from Equation (35). From Table 2 again we have 


d A_ j2 y 2 


(3 - 4r 2 ) 


Hence, we can calculate dQ/dy by means of Equation (36). 


Finally, having calculated dQ/de and dQ/dy we may derive the perturbation 
in co by means of Equation (31). 


32 


ARBITRARY ECCENTRICITY 
TABLE 1 
q = 0, q* = 0 













ARBITRARY ECCENTRICITY 
TABLE 2 
q = 2, q' = - 2 



l-2 e 


C iq (1 +/ 3 2 ) 


-/3 s 

/3 4 (5 + /3 2 ) 

- f3 3 (10 + 5 /3 2 ) 
10 /3 2 (1 + /3 2 ) 

- 5 /3 (1 + 2 /3 2 ) 


- 5 /3 4 + /3 6 

4 /3 3 (5 - /3 2 ) 

5 /3 2 (- 2 + /3 2 ) (3 + /3 2 ) 
20 /3 (1 - /3 4 ) 

- 5 (1 - 2 /3 2 ) (1 + 3 / 3 2 ) 

4 p (l - 5 /3 2 ) 

- 1 + 5/3 2 
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ARBITRARY ECCENTRICITY 

TABLE 3 
q = 2, q' - 0 



c'2 ya 


c' y' (1 - 2 y' 2 ) c y 3 

(1 - 6 y' 2 + 6 y' 4 ) c 2 y : 

-3 c' / (1-2 y' 2 ) c 3 y 
~c' 2 y' 2 c 4 


dA 

Vq 


dy 


4 y 3 

- 2 

c~ 1 y 2 (3 - 4 y 2 ) 

- 1 

2 y (1 - 2 y 2 ) 

0 

c (1 - 4 y 2 ) 

1 

- 4 c 2 y 

2 



— e' 2 
4 

— e' 

2 

(1 - e 2 )~ 3/2 

le' 

2 

le' 2 

4 


c iq (i + /3 2 r 3 

dC . 

iq (1 +/3 2 )" 4 

dy 

-/3 s 

- 5 /3 4 + /3 6 

/? 4 (5 + /3 2 ) 

4 /3 3 (5 - /3 2 ) 

- /3 3 (10 + 5 /3 2 ) 

5 /3 2 (- 2 + /3 2 ) (3 n /3 2 ) 

10 {3 2 (1 + /3 2 ) 

20 /3 (1 - /3 4 ) 

- 5/3 (1 + 2 /3 2 ) 

- 5 (1 - 2 /3 2 ) (1 + 3 /3 2 ) 

1 + 5 /3 2 

4 /S (1 - 5 /3 2 ) 

- P 

- 1 + 5 /3 2 
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ARBITRARY ECCENTRICITY 


TABLE 4 

q = 0, q' = 2 



27(1-2 x 2 ) 

c- 1 (1 - 8 y 2 + 8 y 4 ) 

- 12 x (1-2 -x 2 ) 
c" 1 (1 - 8 y 2 + 8 y 4 ) 


cy(l-2y 2 ) 

1 - 6 y 2 + 6 y 4 
c y (1 - 2 X 2 ) 



-P 3 

- 3 P 2 (1 -P 2 ) 

3 P 2 (1 + P 2 ) 

6 / 3 ( 1 - / 3 4 ) 

- 3 P (1 + 3 P 2 + / 3 4 ) 

- 3 (1 - P 2 ) (1 + 5 P 2 + P 4 ) 

1 + 9 P 2 + 9 / 3 4 + P 6 

12 / 3(1 - / 3 4 ) 

- 3 P (1 + 3 P 2 + P A ) 

- 3 (1 - P 2 ) (1 + 5 / 3 2 + / 3 4 ) 

3 P 2 (1 + P 2 ) 

6/3 (1 - / 3 4 ) 

- P 3 

— 3 /3 2 (1 - / 3 2- ) 


t 
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ARBITRARY ECCENTRICITY 

TABLE 5 
q = 2, q' = 2 


c 2 y 2 


B . - 
i q 


I - _e' 
2 


He ' 2 

2 


dA Vq 


d y 



4 y 3 

- 2 

c -1 y 2 (3 - 4 y 2 ) 

- 1 

2 y (1 - 2 y 2 ) 

0 

c (1 - 4 y 2 ) 

1 

- 4 c 2 y 

2 


C iq (It/? 2 )" 3 


-P s 

/3 4 (5 +/ 3 2 ) 

- P 3 (10 + 5 P 2 ) 
10 P 2 (1 + p 2 ) 

- 5/3 (1 + 2 p 2 ) 
1+5 /3 2 

-P 


(1 +/? 2 )‘ 


-Sp* +P 2 


4/S 3 (5 - /3 2 ) 

5 P 2 (- 2 + £ 2 ) (3 + ^ 2 ) 
20/3 (1 -/3 4 ) 

- 5 (1 - 2 / 3 2 ) (1 + 3 p 2 ) 
4 P (1 - 5 P 2 ) 

-1+5 P 2 
















ARBITRARY ECCENTRICITY 


TABLE 6 
q = 1, q' = 1 



A vq 

dA Vq 

dr 

V 

45 c , 4 ,2 
8 

c 2 y 4 

2 y 3 (2 - 3 y 2 ) 

- 3 

~ c ' 3 (1-3 y' 2 ) / 
4 

c y 3 (2 - 3 y 2 l 

c " 1 y 2 (6 - 23 y 2 + 18 y 4 ) 

- 2 

2c ' 2 (l - 10 / 2 + 15 / 4 ) 
8 

; 2 6 - 20 > 2 + 15 > 4 ) 

2 y (6 - 40 y 2 + 45 y 4 ) 

- 1 

-2c' / (1 - 5 y ' 2 + 5/ 4 ) 

c y (1 - 5 y 2 + 5 y 4 ) 

c ' 1 (1 - 17 y 2 + 45 y 4 - 30 y 6 ) 

0 

2 , '2 (f, - 20 y ' 2 + 15 y’ 4 ) 
8 

c 2 (1 - 10 y 2 + 15 y 4 ) 

- 2 y (11 - 50 > 2 + 45 y 4 ) 

1 

^C' y ' 3 (2-3 y' 2 ) 
4 

c 3 (1 - 3 y 2 ) y 

c (1 - 13 y 2 + 18 y 4 ) 

2 

& 

1*5 

KJ 

c 4 y 2 

2 c 2 y (1 - 3 y 2 ) 

3 


B .V 


i' 


11 e -2 
8 


- 1 


e' 


0 * 


1 + 2 e' 2 


1 


3 e' 


2 


53 e ,j 
8 


3 


C iq (1+/3 2 )- 4 

d ; q < i+ ^ s 

i 


-/3 s 

- /3 4 (5 - 3 /3 2 ) 

- 4 


/3 4 (5 + 3 /3 2 ) 

2 /3 3 (10 - /S 2 - 3 (3 4 ) 

- 3 


- /3 3 (10 + 15 /3 2 + 3 / 3 *) 

-/3 2 (30 + 25 /3 2 - 24 /3 4 - 3 /3 6 ) 

- 2 


(3 2 (10 + 30 /3 2 + 15 /3 4 + /3 6 ) 

2 /3 (10 + 30 /3 2 - 15 /3 4 - 11 /3 6 ) 

- 1 


- 5 /3 (1 + 6 /3 2 + 6 /3* + / 6 ) 

- 5 (1 - /3 4 ) (1 + 11 /3 2 + /3 4 ) 

0 * 


1 + 15 /3 2 + 30 (3* + 10 /3 6 

2 /3 (11 + 15 /3 2 - 30 /3 4 - 10 /3 6 ) 

1 


- /? (3 + 15 /3 2 + 10 /3“) 

- 3 - 2 /3 2 + 25 /? 4 + 30 /S* 

2 


/3 2 (3 + 5 /3 2 ) 

2 /3 (3 + /3 2 - 10 /3 4 ) 

3 


- /3 3 

- /3 2 (3 - 5/S 2 ) 

4 
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ARBITRARY ECCENTRICITY 


TABLK 7 
q =1, q' = - 1 


A 'v q ' 

A Vq 

dA Vq 

dy 

V 

iic' 2 y' 4 
8 

c 2 y 4 

2 y 3 (2 - 3 y 2 ) 


ile' y' 3 (2-3 y' 2 ) 

c y 3 (2 - 3 y 2 ) 

c-> y 2 (6 - 23 y 2 + 18 y 4 ) 


|-y' 2 (6 - 20 y' 2 + IS y' 4 ) 
8 

y 2 (6 - 20 y 2 + 15 y 4 ) 

2 y (6 - 40 y 2 +45 y 4 ) 

- 1 

jC / (1 - S/ ! + S /") 

c y (1 - 5 y 2 + Sy 4 ) 

c* 1 (1 - 17 y 2 + 45 y 4 - 30 y 6 ) 

0 

|-c' 2 (1 - 10 y' 2 + IS y' 4 ) 
8 

c 2 (1 - 10 y 2 + 15 > 4 ) 

- 2 y (11 - 50 y 2 + 45 y 4 ) 

l 

iic-' 3 (1-3 y' 2 ) y 
4 

c 3 (1 - 3 y 2 ) y 

c. (1 - 13 y 2 + 18 y 4 ) 

2 

45 c i4 .2 
8 

c 4 y 2 

2 c 2 y (1 - 3 y 2 ) 

3 

' 

B. , / 
» q 




53 . 2 
8 


- 3 


3 e 


_ 2 


1 + 2 e' 2 


- 1 


e' 


0* 


11 e' 2 

T 


■ 


HKO 

(1 + /3 2 )* 5 

op 

■ 


-/3 s 

- /3 4 (5 - 3 /3 2 ) 

1 


/3 4 (5 + 3/3 2 ) 

2 /3 3 (10 - P 2 - 3 /3 4 ) 

81 


- P 3 (10 + 15 /3 2 + 3 /3 4 ) 

- /3 2 (30 + 25 /2 2 - 24 /3 4 - 3 /? 6 ) 



/S 2 (10 + 30 /3 2 + 15 /3 4 + /3 6 ) 

2 /3 (10 + 30 /? 2 - 15 ,3 4 - 11 /5 6 ) 

- 1 


- 5 /3 (1 + 6 /3 2 + 6,3 4 + £ 6 ) 

- 5 (1 - /3 4 ) (1 + 11 P 2 + /3 4 ) 

0* 


1 + 15 /3 2 + 30 13* + 10 /3 6 

2 ,3 (11 + 15 3 2 - 30 /3 4 - 10 /3 6 ) 

1 


-/3 (3 + 15/5 2 + 10 ^ 4 ) 

- 3 - 24 p 2 + 25 /3 4 + 30 

2 


P 2 (3 + 5 A 2 ) 

2 P (3 + P 2 - 10 /3 4 ) 

3 


- ft 3 

-/ 5 2 (3-5 / 3 2 ) 

4 
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ARBITRARY ECCENTRICITY 



c y 3 ( 2-3 y 2 ) 
y 2 (6 - 20 y 2 + 15 y 4 ) 
c y (1 - 5 y 2 + 5 y 4 ) 
c 2 (1 - 10 y 2 + 15 y 4 ) 
c 3 y (1 - 3 y 2 ) 
c 4 y 2 


2 y 3 (2-3 y 2 ) 

c-> y 2 (6 - 23 y 2 + 18 y 4 ) 

2 y (6 - 40 y 2 + 45 y 4 ) 
c-> (1 - 17 y 2 + 45 y 4 - 30 y 6 ) 
- 2 / (11 - 50 > 2 4 45 y 4 ) 
c (1 - 13 > 2 + 18 ; 4 ) 

2 c 2 y 



1 -6 e ’ 2 
5 e' 


IHZe' 2 



C iq (l + /3 2 )- 4 


-P 5 

P 4 (5 + 3 /3 2 ) 

- ft 3 (10 + 15 P 2 + 3 /3 4 ) 

/3 2 (10 + 30 ft 2 + 15 P 4 + /J 6 ) 

- 5 P (1 + 6 /3 2 + 6 P 4 + (i 6 ) 
1 + 15 /3 2 + 30 p 4 + 10 /3 6 

- p (3 + 15 P 2 + 10 P 4 ) 
p 2 (3 + 5 J 2 ) 

- P 3 


dC iq 

(1 + /3 2 )' 5 

d P 



- /3 4 (5 — 3 P 2 ) 

2 /3 3 (10 - p 2 - 3 /? 4 ) 

- p 2 (30 + 25 « 2 - 24 p 4 - 3 p 6 ) 
2 P (10 + 30 /3 2 - 15 /3 4 - 11 /3 6 ) 

- 5 (1 - P 4 ) (1 + U P 2 +/3 4 ) 

2 P (11 + 15 P 2 - 30 /3 4 - 10 P 6 ) 

- 3 - 24 P 2 + 25 p 4 + 30 p 6 
2 ;i (3 + P 2 - 10 4 ) 

- P 2 (3 - 5 P 2 ) 














ARBITRARY ECCENTRICITY 

TABLE 9 
q = 1, q' = - 3 
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ARBITRARY ECCENTRICITY 


TABLE 10 
q = 3, q' = 1 


A V 

A Vq 

d A y q 

d / 

Bi 

lie ' 4 y' 2 
8 

y 6 

6 y 5 


4 c ' 3 / (i - 3 y 2 ) 

4 

c ; s 

c " 1 y 4 (5 - 6 y 2 ) 


lie -' 3 ( 1-10 y' 2 + 15 y M ) 
8 

c J y 4 

2 / 3 (2 - 3 y 2 ) 

- 1 

(1-3 ' 2 + 5 /*) 

c 3 y 3 

3 c / 2 (1 - 2 y 2 ) 

0 

4 >'* (6 - 20 >' 2 + IS y' 4 ) 

O 

C 4 y s 

2 c 2 / (1 - 3 y 2 ) 

1 

-~c y ' 3 (2-3 y' 2 ) 
4 

c 5 y 

c 3 (1 - 6 y 2 ) 


15 ,2 . 1 4 

C ^ 

8 

y 

- 6 c 4 > 



B * • 
» q 


| 


11 e ,2 
8 
e' 


- 1 
0 * 


1 + 2 e ' 2 


1 


3 e' 


2 


53 e , 2 
8 


3 



dC. , , 

iq (i + P 2 y s 

dp 

i 


-p 7 

- P b (7 - P 2 ) 

- 4 


P 6 (7 + P 2 ) 

6 /3 s (7 - P 2 ) 

- 3 


-7/35 ( 3 +/ 8 2 ) 

- 7 p* (5 + P 2 ) (3 - P 2 ) 

- 2 


7 /S 4 (5 + 3 /3 2 ) 

14 P 3 (10 - P 2 - 3 /3 4 ) 

- 1 


- 35 /3 3 (1 + /3 2 ) 

- 105 A ? 2 (1 - P 4 ) 

0 * 


1 ft 2 (3 + Sp 2 ) 

14 P (3 + , 6 2 - 10 /3 4 ) 

1 


-7/6(1 +3 P 2 ) 

- 7 (1 - 3 /3 2 ) (1 + 5 p 2 ) 

2 


1 + 7 p 2 

6/3(1 - 7 P 2 ) 

3 


- P 

- (1 - 7 /3 2 ) 

4 


t 
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ARBITRARY ECCENTRICITY 


TABLE 11 
q = 3, q' si.- 1 
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ARBITRARY ECCENTRICITY 


TABLE 12 



q = 3, q 

' = 3 


A vq* 

A vq 

d A V<1 
d7 

V 

ic * 6 

8 

y‘ 

6 7 s 

- 3 

- x 1 '* >■ 

4 

c z s 

c ' 1 > 4 (5 - 6 7 2 ) 

- 2 

75 c 1 4 yj 
8 

e 2 7 4 

2 / 3 (2-3 y 2 ) 

- 1 

-He' 3 /* 
2 

c 3 7 3 

3 c y 2 (1 - 2 x 2 ) 

0 

75 c <2 y« 
8 

c 4 7 2 

2 e 2 7 (1 - 3 7 2 ) 

1 

15 c < l 5 
4 

c 5 > 

c 3 (1 - 6 7 2 ) 

2 

f '’ 1 

c« 

- 6 e 4 > 

3 


V,' 


i * 


ie ' 2 

8 

- e' 


1 

2 


1 - 6 e ’ 2 


3 


5 e' 


4 


127 ..a 


5 


8 




C iq (1 + /32) ’ 4 

dC lq (1 +/>*)-* 

d/3 

i 


-/3 7 

-/ 3 6 (7 -/ 3 2 ) 

- 4 


/3 6 (7 + /3 2 ) 

6 /3 s (7 -/3 2 ) 

- 3 


- 7 /3 5 (3 + /3 2 ) 

- 7 /3 4 (5 +/3 2 ) (3 - / 3 2 ) 

- 2 


7 /3 4 (5 + 3/3 2 ) 

14 /3 3 (10 - , 5 2 - 3 / 3 4 ) 

- 1 


- 35 /3 3 (1 + /3 2 ) 

- 105 /3 2 (1 -/ 3 4 ) 

0 


7 /3 2 (3 + 5 /3 2 ) 

14 /3 (3 +/3 2 - 10 /3 4 ) 

1 


- 7 /3 (1 + 3 /3 2 ) 

- 7 (1 - 3 / 3 2 ) (1 + 5 / 3 2 ) 

2 


1 + 7 / 3 2 

6 /3 (1 - 7 / 3 2 ) 

3 


-/* 

- (1 - 7 /3 2 ) 

4 
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ARBITRARY ECCENTRICITY 









TABLE 14 (Continued) 


« 


^2 _ i' 2 (i - 2) cos [(i - 2) E + i* / t* + 0] i’ 3 4 cos (i E + i* + <j>) 

4[(i-l) + i'm][(i-2) + i'm] 2[(i-l) + i'm][(i+l) + i'm] 


<# ” 

i' 2 (i + 2) cos [(i +2) E + i' l' + <£] 
4[(i+l) + i'm] [(i + 2) + i' m] 


^2 qi' 2 cos [(i-2)E+i , 't , t^] qi' 2 cos(iE + i' -£'+<£) 

3ei ' 4 [ (i -1) + i'm] [(i -2) + i'm] 2 [(i - 1) + i' m] [(i + 1) + i' m] 


q i' 2 cos [(i + 2) E + i' i' + <p] 
4 [(i +1) + i'm] [(i + 2) + i'm] 


_ ^i' 2 cos [(i -2) E + i' 'fc' + <£] + v i' 2 cos (i E + i'l' + 0) 
dd 4 [(i -1) + i'm] [ (i -2) + i' m] 2 [(i - 1) + i ' m] [(i + 1) + i' m] 


i^i' 2 cos [(i +2) E+ i' V + <p] 
4[(i+l) + i'm] [(i + 2) + i' m] 


^ 3 i' (i + 1) me sin [(i + 1) E + i' / t' + <£] 

^ a 4 a [(i + 1) + i' m] 2 


3 i* (i - 1) me sin [(i— 1)E + i' l' +0] 

4 a [(i -1) + i'm] 2 
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TABLE 14 (Continued) 

a [m 2 e 2 Q'] _ 3 j<2 m 2 e 2 [ 2 (i -1) (i -2) + i*(2 i -3) m] sin [(i - 2) E + i *1* + <f>] r 

8a [(i-1) + i'm] 2 [ (i — 2) + i'm] 2 

" * 


3 i' 2 m 2 e 2 [i 2 - 1 + m i i'] sin (i E + i' •£' + 0) 
2a [(i -1) + i'm] 2 [(i + 1) + i'm] 2 


3i* 2 m 2 e 2 [2 (i + 1) (i + 2) + i * (2 i + 3) m] sin [(i + 2)E + i*4‘ +<£] 
8a [(i + 1) + i'm] 2 [(i-2) + i'm] 2 
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